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Abstract 

A short discussion of theoretical results on cosmic ray first-order Fermi acceleration at relativistic shock waves is 
presented. We point out that the recent results by Niemiec with collaborators change the knowledge about these 
processes in a substantial way. In particular one can not expect such shocks to form particle distributions extending 
to very high energies. Instead, distributions with the shock compressed injected component followed by a more or less 
extended high energy tail are usually created. Increasing the shock Lorentz factor leads to steepening and decreasing 
of the energetic tail. Also, even if a given section of the spectrum preserves the power-law form, the fitted spectral 
index may be larger or smaller than the claimed 'universal index' a « 2.2 . 

A reported simple check of real shapes of electron spectra in the Cyg A hot spots provides results clearly deviating 
from the standard expectations for such shocks met in the literature. The spectrum consist of the very fiat low energy 
part (a « 1.5), up to electron energies ~ 1 GeV, and much steeper part [a > 3) at higher energies. We conclude with 
presentation of a short qualitative discussion of the Fermi second-order processes acting in relativistic plasmas. We 
suggest that such processes can be the main accelerating agent for very high energy particles. In particular its can 
accelerate electrons to energies in the range of 1 - 10'^ TeV in relativistic jets, shocks and radio-source lobes. 
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1. Introduction 

Relativistic plasma flows are observed in a number 
of astrophysical objects, ranging from a mildly rel- 
ativistic jets of the sources like SS433, through the- 
Lorentz-factor-of-a-few jets in AGNs and galactic 
'micro-quasars', up to highly ultra-relativistic out- 
flows in sources of gamma ray bursts or pulsar winds. 
As nearly all such objects are eflicient emitters of 
nonthermal radiation, what requires existence of en- 
ergetic particles, our attempts to understand the 
processes generating cosmic rays are essential for 
understanding the fascinating phenomena observed. 
Below I will discuss the work carried out in order 
to understand the cosmic ray acceleration processes 
acting at relativistic shocks and within highly tur- 



bulent regions accompanying such shocks and shear 
layers. I will not include here the interesting work 
involving collisionless shocks modelling with par- 
ticle in cell simulations. This approach uses quite 
different modelling method in comparison to the 
other work discussed here, relating in most cases to 
the characteristic energy range of shock compressed 
thermal plasma particles. 

The present paper is a modified and updated ver- 
sion of some of my previous reviews of the subject. 
Also, essentially the same slightly shortened text is 
provided as my contribution to the HEPRO Work- 
shop in Dublin (September 2007). 

Below we will append the index '1' ('2') to quan- 
tities measured in the plasma rest frame upstream 
(downstream) of the shock. 
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2. Particle diffusive acceleration at 
non-relativistic shock waves 

Processes of the first-order particle acceleration 

at non-rclativistic shock waves were widely dis- 
cussed by a number of authors, let me note still 
actual reviews by Drury (1983) and Blandford & 
Eichlcr (1987). The most interesting physical fea- 
ture of the first-order shock acceleration at the 
non-relativistic shock wave is independence of the 
test-particle stationary particle energy spectrum 
from the background conditions near the shock, in- 
cluding the mean magnetic field configuration and 
the spectrum of MHD turbulence. The reason is 
a nearly-isotropic form of the particle momentum 
distribution at the shock. If efficient scattering oc- 
curs near the shock, this condition also holds for 
oblique shocks with the shock velocity along the 
(upstream) magnetic field Ub,i = U\/ cos'J'i <C v 
{^i - the upstream magnetic field inclination to the 
shock normal). Then, the particle density is contin- 
uous across the shock and the spectral index for the 
phase-space distribution function, a, is given exclu- 
sively in terms of a single parameter - the shock 
compression ratio R: 



Because of the isotropic form of the particle distri- 
bution function, the spatial diffusion equation has 
become a widely used mathematical tool for describ- 
ing particle transport and acceleration processes in 
non-relativistic flows. 

In real astronomical objects some non-linear and 
time dependent effects, or occurring of additional 
energy losses and gains can make the physics of the 
acceleration more complicated, creating, e.g., non- 
power- low and/or non-stationary particle distribu- 
tions. 

3. Cosmic ray acceleration at relativistic 
shock waves 

Below I describe work done on mildly- and ultra- 
relativistic shock acceleration including important 
recent results of Niemiec et al. With these last re- 
sults many previous ones occurred to be of histori- 
cal value only, reflecting specific individual features 
of the acceleration processes. Basing on these older 
works one can understand the recent simulations re- 
sults in a relatively straightforward way. Attempt- 



ing to give an overview of the full field I base this 

presentation on my own and my collaborators work, 
which seems to me to present a consistent way of de- 
velopment and reflects my approach to understand- 
ing acceleration processes in relativistic shocks. 

3.1. History: acceleration at mildly relativistic 
shocks 

In cases of the shock velocity reaching values com- 
parable to the light velocity, the particle distribu- 
tion at the shock becomes anisotropic. This simple 
fact complicates to a great extent both the physical 
picture and the mathematical description of particle 
acceleration. A first attempt to consider the acceler- 
ation process at the relativistic shock was presented 
in 1981 by Peacock, and a consistent theory was 
proposed later by Kirk & Schneider (1987a). Those 
authors considered stationary solutions of the rel- 
ativistic Fokker-Planck equation for particle pitch- 
angle diffusion in the parallel shock wave. In the 
situation with the gyro-phase averaged distribution 
/(p, /X, z), which depends only on the unique spatial 
co-ordinate z along the shock velocity, and with fi 
being the pitch-angle cosine, the equation takes a 
form: 

T{U + v^l)^^=C{f) + S , (3.1) 

where F = — is the flow Lorentz factor, 

C{f) is the collision operator and S is the source 
function. In the presented approach, the spatial co- 
ordinates are measured in the shock rest frame, while 
the particle momentum co-ordinates and the col- 
lision operator are given in the respective plasma 
rest frame. For the applied pitch-angle diffusion op- 
erator, C = d/dii{D^^df /djj), the authors gener- 
alized the diffusive approach to higher order terms 
in particle distribution anisotropy and constructed 
general solutions at both sides of the shock which 
involved solutions of the eigenvalue problem. By 
matching two solutions at the shock, the spectral in- 
dex of the resulting power-law particle distribution 
can be found by taking into account a sufBciently 
large number of eigenfunctions. The same procedure 
yields the particle angular distribution and the spa- 
tial density distribution. The low-order truncation 
in this approach corresponds to the standard diffu- 
sion approximation and to a somewhat more general 
method described by Peacock. 
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An application of this approach to more reahstic 
conditions - but still for parallel shocks - was pre- 
sented by Heavens & Drury (1988), who investigated 
the fluid dynamics of relativistic shocks (cf. also El- 
lison & Reynolds 1991) and used the results to cal- 
culate spectral indices for accelerated particles. 
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Fig. 1. Spectral indices a of particles accelerated at oblique 
shocks versus shock velocity projected at the mean magnetic 
field, Ub,i (Ostrowski 1991a). The results are presented for 
the shock compression R = 4. On the right the respective 

synchrotron spectral index 7 is given. The shock velocities 
Ui are given near the respective curves taken from Kirk 
& Heavens (1989). The points were taken from simulations 
computing explicitly the details of particle-shock interactions 
(Ostrowski 1991a). 

A substantial progress in understanding the accel- 
eration process in the presence of highly anisotropic 
particle distributions is due to the work of Kirk 
& Heavens (1989; see also Ostrowski 1991a and 
Ballard & Heavens 1991), who considered particle 
acceleration at subluminal {Ub,i < c) relativistic 
shocks with oblique magnetic fields. They assumed 
the magnetic momentum conservation, p\/B = 
const, at particle interaction with the shock and ap- 
plied the Fokker-Planck equation discussed above 
to describe particle transport along the field lines 
outside the shock, while excluding the possibility of 
cross-field diflfusion. In the cases when Ubs reached 
relativistic values, they derived very flat energy 
spectra with 7 « at Ub,i ~ 1 (Fig- !)• In such 
conditions, the particle density in front of the shock 
can substantially - even by a few orders of mag- 
nitude - exceed the downstream density (see the 
curve denoted '-8.9' at Fig. 2). Creating flat spec- 
tra and great density contrasts is possible due to 
the effective repeating reflections of anisotropically 
distributed upstream particles from the region of 
compressed magnetic fleld downstream of the shock. 



As stressed by Begelman & Kirk (1990), in rel- 
ativistic shocks one often flnds superluminal con- 
ditions with Ub,i > c, where the above presented 
approach is no longer valid. Then, it is not possi- 
ble to reflect upstream particles from the shock and 
to transmit downstream particles into the upstream 
region. In effect, only a single transmission of up- 
stream (or shock injected) particles re-shapes the 
original distribution by shifting particle energies to 
larger values, with super-adiabatic efficiency due to 
anisotropic particle distribution at the transmission. 
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Fig. 2. Energetic particle density profiles across the rel- 
ativistic shock with the oblique magnetic field (Ostrowski 
1991b). The shock with Ui = 0.5, R = 5.11 and Vi = 55° 
is considered. The curves for growing to the top perturba- 
tion amplitudes are characterized with the value logKx/K|| 
{kj_/k^^ is a ratio of the cross-field to the parallel diffusion 
coefficients) given near each curve. The data arc vertically 
shifted for picture clarity. The value Xmax is the distance 
from the shock at which the upstream particle density de- 
creases to 10~^ part of the shock value. 

3.2. History: acceleration in the presence of large 
amplitude magnetic field perturbations 

As the approaches proposed by Kirk & Schnei- 
der (1987a) and Kirk & Heavens (1989), and the 
derivations of Begelman & Kirk (1990) are valid 
only in cases of weakly perturbed magnetic fields, 
for larger amplitude MHD perturbations numerical 
methods have to be used. A series of such simu- 
lations were performed by numerous authors (e.g. 
Kirk & Schneider 1987b; ElHson et al. 1990; Os- 
trowski 1991a, 1993; Ballard & Heavens 1992, Naito 
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& Takahara 1995, Bednarz & Ostrowski 1996). Dif- 
ferent approaches applied, including the ones involv- 
ing particle momentum pitch-angle diffusion or in- 
tegrating particle trajectories in analytic structures 
of the perturbed magnetic fieldQ, allowed for con- 
sidering a wide range of background conditions near 
the shock. The results obtained by different authors 
can be summarized at the figure (Fig. 3) taken from 
Bednarz & Ostrowski (1996). One should note, that 
essentially all derivations by the above authors were 
performed with assuming scale-free conditions for 
the acceleration process, resulting in power-law dis- 
tributions of the accelerated particles. 



\ 


Ui = 0.5 


2 


5.8° 




45.6° 

I 






\ \ , 


^~~~~~~-~_______72.5° 




60° " ~* 



^ 3.0 3.5 4.0 4.5 5.0 5.5 

a 



Fig. 3. Relation of the particle acceleration time scale Tacc 
versus the particle spectral index a at different magnetic 
field inclinations ipi given near the respective curves (Bed- 
narz & Ostrowski 1996). The minimum value of the model 
parameter occurs at the encircled point of each curve 

and the wave amplitude monotonously increases along each 
curve up to SB ~ B, where all curves converge in a. The 
curve for = 60° {Ub,i = 1) separate the sub- and super- 
luminal shock results. We do not discuss here the presented 
acceleration times. 

At the figure (Fig. 3) one can find very flat spec- 
tra for oblique subluminal shocks if the perturba- 
tion amplitudes are small. Contrary to that genera- 
tion of the power-law spectrum is possible in the su- 
perluminal shocks only in the presence of large am- 
plitude turbulence. Then, in contrast to the sublu- 
minal shocks, spectra are extremely steep for small 
values of SB (not presented at the figure) and a 
monotonously decreases with increasing magnetic 
field perturbations. A new feature is observed in 



^ Let us note that application by some authors of point-like 
largo-angle scattering models in relativistic shocks does not 
provide a viable physical representation of the scattering at 
MHD waves (Bednarz & Ostrowski 1996). 



oblique shocks of the spectral index a changing with 
SB in a non-monotonic way. 

3.3. History: Energy spectra of cosmic rays 
accelerated at ultra-relativistic shocks 

The main difficulty in modelling the acceleration 
processes at shocks with large Lorentz factors F is 
the fact that the involved particle distributions are 
extremely anisotropic in the shock, with the parti- 
cle angular distribution opening angles ^ in the 
upstream plasma rest frame. In the simulations of 
Bednarz & Ostrowski (1998) a Monte Carlo method 
involving small amplitude pitch-angle scattering was 
applied for particle transport near the shocks with F 
in the range 3 - 243. The simulations revealed an un- 
expected result, showing convergence, for F — > oo, 
of the resulting power-law distributions to the 'uni- 
versal' one with the spectral index a « 2.2 (Fig. 4). 
Essentially the same result was derived with dif- 
ferent methods by many other authors (Gallant & 
Achterberg 1999; Kirk et al. 2000, Achterberg et al. 
2001, Lemoione et al. 2003, Keshet & Waxman 2005, 
Lemoine & Revenu 2006, Morhno et al. 2007, et al.), 
what could suggest that in the ultrarelativistic limit 
the acceleration process becomes again simple, gen- 
erating cosmic ray spectra essentially independent 
of the considered background conditions. One can 
find examples in the literature that some authors 
extend this claim for all relativistic shocks. 









'' 30°'- 


90° 




X= -3.44 




20°/' 


', 60°^'. 








10° 




















, - ' 0° 


& 


: ;/|;/ 


; > i i i 


3 


9 


27 

r 


1 

81 


1 

243 



Fig. 4. The simulated spectral indices a [a = a — 2) versus 
the shock Lorentz factor P (Bednarz & Ostrowski 1998). 
Results for a given are joined with dashed lines; the 
respective value of ipi is given near each curve. Increasing 
the turbulence amplitude in a not presented here series of 
simulations led to shifting the resulting curves toward the 
parallel shock, = 0°, results. 
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Fig. 5. Particle spectra for an oblique mildly relativistic shock: the shock velocity ui = 0.5 c, the mean magnetic field inclination 
^1 = 45° and the wave power spectrum are indicated in the respective panels (from Niemiec & Ostrowski 2004). Values of 
magnetic turbulence amplitude, 5B/B, and the indices fitted to the power-law sections of spectra (in parentheses) are given 
near each result. 



Ostrowski & Bednarz (2002) reconsidered all the 
above approaches to derive particle spectra at rela- 
tivistic shocks and 'discovered' that the conditions 
producing the universal spectral index were in some 
way equivalent to assuming highly turbulent condi- 
tions near the shock. Additionally, all these mod- 
els did not introduce any physical scale and thus 
forced the power-law shape of the resulting spec- 
trum. Do such conditions and the resulting charac- 
teristic spectra really exist in astrophysical situa- 
tions ? 



3.4. Toward a realistic description oj the relativistic 
shock acceleration 

Studies of - as far as possible - realistic conditions 
for the relativistic shock acceleration were presented 
in a series of papers by Niemiec et al. (Niemiec & 
Ostrowski 2004, 2006, Niemiec et al. 2006; see also 
Lemoine Pelletier 2006). We assumed 3-D static 
magnetic field perturbations upstream of the shock 
by imposing a large number of sinusoidal waves 
with different power spectra, F{k). The flat spec- 
trum, F{k) cx , or the Kolmogorov spectrum, 
F{k) cx k~^^^, were considered in the wide wave- 
vector range (kmim ^ma£c); with k„iax/kmin — 10^ 

The downstream magnetic field structures were 
computed by respective compression of the up- 
stream field at the shock. In the last of above papers 



an additional component of large amplitude short- 
wave MHD turbulence was assumed to be produce 
at the shock. The accelerated particle spectra were 
derived using the Monte Carlo simulations for a 
wide range of shock Lorentz factors - between 2 
and 30 - and a selection of the mean magnetic field 
configurations and perturbation amplitudes. 
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Fig. 6. The particle spectra derived for superluminal rel- 
ativistic shocks with Lorentz factors 71 = 5, 10 and 30 
(Niemiec & Ostrowski 2006a). 

Let us take a look at a few characteristic results 
of these modelings. At Fig. 5 results derived for sub- 
luminal oblique mildly relativistic shocks are pre- 



5 



O 



O 



7 I 
6 

5 

4 - 
3 - 
2 - 
1 - 
- 
-1 - 



a -3. 6 



7, = 30 C'^) 

-^1 = 0° 
F(k)~k"' - 



\ 



6B/Bo_,=0.05 



a 9-5 

ft 



27T/k„ 



a'VXA 0.1 
\ A 

A 

X X 3.0 
□ 1.0 
+ + 0.3 



A 



a-3.5 



7, = 30 C^) 

V,=o° 

F(k)~k-^/^- 



(5B/Bo,,=0.05 



A 0.1 

A 



'•ta 



27T/k„ 







2 5 4 5 6 7 1 

log(E/Eo) 



2 5 4 5 6 7 



Fig. 7. Particle spectra derived for parallel relativistic shocks with Lorentz factor 71 = 30 (Nicmicc & Ostrowski 2006). For 
description see Fig. 5 and the original paper. 



sented. One may note that introducing the energy 
scales to the modelHng, in our units 2t: /kmax and 
2T:/kmim lead to deviations of the resulting spectra 
from the power-law shape. The spectra, as expected 
from the previous discussion, are very flat for small 
amplitude turbulence, but steepen at larger ampli- 
tudes. An interesting feature is seen, that at very 
high particle energies, above the resonance range 
{E > 27r/fcmi„) the 'short' waves weaker influence 
particle trajectories leading to the hard energy tails 
before the cut-offs imposed by the modelling. 

When we considered oblique superluminal shocks 
the spectra consisted of the shock compressed in- 
jected component and the limited high energy tail 
until a cut-off well within the 'resonance range' 
{2Tr/kmax, 27r/fcmin). As illustrated at Fig. 6 the 
tails fast diminish with the growing shock Lorentz 
factor, leaving for large 71 the 'compressed compo- 
nent' only. 

There are a few general observations for the first- 
order Fermi acceleration processes from these series 
of models. For particles in a low energy range of the 
resonance energies for the considered field pertur- 
bations the acceleration processes proceed in an en- 
semble of different oblique shocks, where each local 
mean magnetic field structure is formed as a super- 
position of the mean magnetic field i?o,i and long 
wave field perturbations. Thus, there can occur sig- 
nificant differences between spectra generated in the 
presence of flat and steep turbulence power spec- 
tra, and the spectral indices significantly vary with 



the perturbations amplitudes. In parallel shocks the 
long wave perturbations introduce the acceleration 
effects observed with oblique magnetic fields (cf . also 
Ostrowski 1993). Thus in the ultrarelativistic paral- 
lel shocks propagating in highly turbulent medium 
these effects can lead to formation of particle dis- 
tributions with cut-offs at relatively low energies, 
like in shocks with perpendicular field configura- 
tions (Fig. 7). If the shock wave generates a large 
amplitude short-wave turbulence downstream of the 
shock, the acceleration process can form a more ex- 
tended power-law tail, but at higher particle energies 
the mean magnetic field or the long wave magnetic 
field structures start to dominate in shaping par- 
ticle trajectories and thus the acceleration process, 
leading to results like the ones described above. In 
any studied case we were not able to create the scale 
free conditions in the acceleration process, leading 
to the wide range power law distribution of acceler- 
ated particles. It was possible only in limited energy 
ranges and the forms of the spectra depended usu- 
ally a lot on the considered background conditions. 

3.5. Observational constraints on the shock 
acceleration from Cyg A hotspots 

Constraints for the above theoretical derivations 
can be provided by precise observations of energetic 
particle emission from objects harbouring the rela- 
tivistic shocks. Such study performed for hotspots 



6 



of the Cygnus A radiosource (Stawarz et al. 2007) 
reveals significant deviations of the derived electron 
spectra from the 'standard' shock spectra. The re- 
sulting spectral energy distribution of the hotspot 
D is presented at Fig. 8, showing both the extended 
synchrotron component and the inverse-compton 
(IC) one modelled for optical and X-ray measure- 
ments. Additionally, the low Spitzer IR points 
provide additional important constraint for the IC 
spectral component. In derivation of the rclativistic 
electrons distribution these measurements allow for 
excluding a possibility of substantial absorption in 
the low frequency synchrotron spectrum and thus 
require very flat distribution of low energy elec- 
trons. Thus the intrinsic electron spectrum (Fig. 9) 
is composed of a very flat low energy sector, with 
the energy spectral index s ~ 1.5, followed above 
a break at i? ~ 1 GeV with the steep (s > 3) 
high energy sector. A possible interpretation of the 
spectrum considered by us was to assume a mildly- 
relativistic shock acceleration in the jet dominated 
by the protons bulk kinetic energy. The protons are 
expected establish the characteristic break energy 
scale, Ebr ^ 1 GeV, at the shock transition layer. 
The electrons (or pairs) below this energy are ex- 
pected to be accelerated within this layer due to 
electron-proton collective interactions. Above Ebr, 
either the first order Fermi process acting at the 
shock creates a steep spectrum, or the accelera- 
tion process proceeds downstream of the shock in 
the second-order Fermi process. The existing the- 
oretical models do not allow to reject any of these 
alternatives. 

4. Energetic particle acceleration in shear 
layers and regions of relativistic turbulence 

Acceleration processes acting, e.g., in AGN cen- 
tral sources and in shocks formed in large scale 
jets are not always able to explain the observed 
high energy electrons radiating away from the cen- 
tre/shock. Among a few proposed possibilities ex- 
plaining these data the relatively natural but still 
unexplored is the one involving particle acceleration 
within a shear layer transition at the interface be- 
tween the jet and the surrounding medium. To date 
knowledge of physical conditions within such layers 
is very limited and only rough estimates for the con- 
sidered acceleration processes are possible. Within 
the subsonic turbulent layer with a non-vanishing 
small velocity shear the ordinary second-order Fermi 
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Fig. 8. Spectral energy distribution of emission from the hot 
spot D of Cyg A (Stawarz et al. 2007). One can clearly see 
both the synchrotron and the IC components. The Spitzer 
points in infrared and the optical point allow to exclude 
possibility of forming the measured flat low-frequency syn- 
chrotron component (dotted lines in synchrotron and IC 
spectral ranges) due to some self-absorption processes. 
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Fig. 9. Relativistic electron spectra in Cyg A hot-spots 
(Stawarz et al. 2007; here 7 is an electron Lorentz factor). 
Different spectral indices of the low energy and the high 
energy parts are expected to be intrinsic to the acceleration 
process, not the effect of the distribution 'aging' downstream 
of the shock. 

acceleration, as well as the process of 'viscous' par- 
ticle acceleration (cf. the review by Berezhko (1990) 
of the work done in early 80-th; Earl et al. 1988, 
Ostrowski 1990, 1998, 2000, Stawarz & Ostrowski 
2002) can take place. A mean particle energy gain 
per scattering in the later process scales as 

14^.(1^)' , 

where < AC/ > is the mean velocity difference be- 
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tween the 'successive scattering acts'. It is propor- 
tional to the mean free path normal to the shear 
layer, A„, times the mean flow velocity gradient in 
this direction V„ • U. With d denoting the shear layer 
thickness this gradient can be estimated as | V„-C/| « 
U/d. Because the acceleration rate in the Fermi II 
process is cx {V/c)^ {V « Va is the wave velocity, 
Va - the Alfven velocity), the relative importance 
of both processes is given by a factor 

my ■ 

The relative efficiency of the viscous acceleration 
grows with A„ and in the formal limit of A„ « d and 
1/ <C c outside the equation (4.2) validity range it 
dominates over the Fermi acceleration to a large ex- 
tent. Because accelerated particles can escape from 
the accelerating layer only due to a relatively inef- 
ficient radial diffusion, the formed particle spectra 
are expected be very flat up to the high energy cut- 
off, but the exact form of the spectrum depends on 
several unknown physical parameters of the bound- 
ary layer (Ostrowski 1998, 2000). 

For turbulent relativistic plasmas the second- 
order Fermi acceleration can in principle dominate 
over the viscous process at all particle energies. In 
the case of electrons the upper energy scale for the 
accelerated particles is provided by the radiation 
losses. A simple exercises with the above estimated 
acceleration scales and the synchrotron radiation 
loss scale yields - for the sources like small and large 
scale jets, radio hot spots or lobes - the highest elec- 
tron energies between 1 TeV and 10^ TeV, for the 
respective acceleration time scales ~ 10^ s in gauss 
fields and up to ~ 10^ yrs in fields (depending 
on the considered object). The much higher energy 
limits for protons are usually determined by the es- 
cape boundary conditions, not the radiative losses. 

5. Final remarks 

A recent study of the first order Fermi accelera- 
tion processes at relativistic shocks, taking into ac- 
count realistic assumptions about the physical con- 
ditions near the shock, reveals a few unexpected con- 
clusions. The modelling of the acceleration process 
in mostly perpendicular (for relativistic velocities) 
shocks yields spectra consisting of the compressed 
injected part appended with a limited high energy 
tail. For given upstream conditions increasing the 
shock Lorentz factor F leads to steepening of such 



energetic tails, providing essentially the compressed 

injected component at large F. Another unexpected 
feature is observed dependence of the spectrum incli- 
nation on the turbulence amplitude also for the par- 
allel shock waves and formation of cut-offs at such 
shocks for large F. Essentially no conditions stud- 
ied by Niemiec et al. resulted in formation of the 
wide-range power-law particle distribution with the 
universal spectral index 2.2 . 

It can be of interest, with the presently published 
AUGER results (The Pierre Auger Collaboration 
2007), that the modelling presented above seems to 
exclude the first-order Fermi acceleration at rela- 
tivistic shocks as possible sources for highest energy 
particles registered in this experiment. 

In the same time the second-order Fermi pro- 
cesses acting in turbulent relativistic plasmas are 
expected to play significant role in cosmic ray accel- 
eration. Thus, even facing substantial mathematical 
and physical difhculties, its deserve a detailed study. 
One may note that the acceleration processes accom- 
panied the magnetic field reconnection processes are 
analogous to the second-order Fermi acceleration. 
Such processes always accompany the large ampli- 
tude MHD turbulence and generate turbulence. A 
few simple attempts to consider these processes were 
recently presented (e.g. Virtanen & Vainio 2005). 

When considering the relativistic shock acceler- 
ation one should also note interesting approaches 
by Derishev et al. (2003) and Stern (2003), outside 
the classical Fermi scheme. They consider the ac- 
celeration processes in highly relativistic shocks or 
jet shear layers (Stern & Poutanen 2006) involving 
particle-particle or particle-photon interactions at 
both sides of the shock. 
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